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ABSTRACT 

The origin of the high-energy component in spectral energy distributions (SED) of blazars 
is still a bit of a mystery. While BL Lac objects can be rather successfully modeled within 
the one-zone synchrotron self-Compton (SSC) scenario, the SED of low peaked Flat Spectrum 
Radio Quasars (FSRQ) is more difficult to reproduce. Their high-energy component needs the 
abundance of strong external photon sources, giving rise to stronger cooling via the inverse 
Compton channel, and thusly to a powerful component in the SED. Recently, we were able to show 
that such a powerful inverse Compton component can also be achieved within the SSC framework. 
This, however, is only possible if the electrons cool by SSC, which results in a non-linear process, 
since the cooling depends on an energy integral over the electrons. In this paper we aim to 
compare the non-linear SSC framework with the external Compton (EC) output by calculating 
analytically the external Compton component with the underlying electron distribution being 
either linearly or non-linearly cooled. Due to the additional linear cooling of the electrons with 
the external photons, higher number densities of electrons are required to achieve non-linear 
cooling, resulting in more powerful inverse Compton components. If the electrons initially cool 
non-linearly, the resulting SED can exhibit a dominating SSC over the EC component. However, 
this dominance depends strongly on the input parameters. We conclude that with the correct 
time-dependent treatment the SSC component should be taken into account to model blazar 
flares. 

Subject headings: radiation mechanisms: non-thermal - BL Lacertae objects: general - gamma-rays: 
theory 



1. Introduction 

Combined as blazars, flat spectrum radio 
quasars (FSRQ) and BL Lacertae objects (BL 
Lacs) arc the most violent subgroup of active 
galactic nuclei (AGN) from the earth's point of 
view in the accepted uniflcation scheme (Urry & 
Padovanni 1995). The broadband spectral energy 
distribution of blazars is dominated by two broad 
non-thermal components. The low-energetic one, 
peaking usually between the infrared and the X- 
ray parts, is attributed to synchrotron radiation 
of highly relativistic electrons, while the process 
behind the high-energy component, peaking in the 



7-rays, is a matter of ongoing discussions. 

Albeit the possibility of a hadronic origin (e.g. 
Mannheim 1993) for the high-energetic compo- 
nent, most authors favor a leptonic origin of the 
7-radiation (for recent reviews see Bottcher 2007, 
2012). If highly relativistic electrons (and also 
positrons) interact with an ambient photon field, 
the photons can be inverse Compton scattered to 
very high energies. Such photon fields can be the 
synchrotron photons of the same population of 
electrons (the so-called synchrotron self-Compton 
effect (SSC), Jones, O'Defl & Stein 1974), or pho- 
ton sources external to the jet (so-called external 
Compton models (EC)), like photons directly from 
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the accretion disk surrounding the black hole in 
the center of the active galaxy (Dermer & Schlick- 
eiser 1993), from the broad line regions (Sikora et 
al. 1994) or the dusty torus (Blazejowski et al. 
2000, Arbeiter et al. 2002). 

Especially the SEDs of FSRQs are dominated 
by the inverse Compton component, i.e. most of 
the luminosity of this type of blazars is emitted 
in 7-rays (e.g. Hayashida et al. 2012 for 3C 279, 
or Vercellone et al. 2011 for 3C 454.3). Since one 
can measure several thermal emission components 
in FSRQs, the EC process seems to be a natu- 
ral choice to model the SED. Being able to detect 
unbeamed thermal emission, albeit the strongly 
boosted non-thermal radiation of the jet, proves 
the abundance of strong external photon fields. 
This in turn provides lots of seed photons for the 
electrons, which cool rather strongly by this pro- 
cess, giving rise to a powerful inverse Compton 
component. 

On the other hand, BL Lacs emit most of their 
power in the synchrotron component having at 
most comparable inverse Compton fluxes. Es- 
pecially high-frequency peaked BL Lacs exhibit 
a much reduces 7-ray flux compared to the syn- 
chrotron flux. This points towards a main cooling 
by the synchrotron channel, and the SSC process 
is successfully used to model BL Lacs (e.g. Acciari 
et al. 2011 for lES 2344+514, or Abramowski et 
al. 2012 for IRXS J101015.9-311909). 

However, such a strict division of the two blazar 
types has been called into question, recently. Chen 
et al. (2012) performed a numerical analysis of 
the multiwavelength variability of the FSRQ PKS 
1510-089. Using a time-dependent code they were 
unable to find a clear preference of the EC over the 
SSC process, and state that the SSC process might 
even be preferable, since it matches the X-ray part 
of the SED far better than the EC process. 

The advantage of the time-dependent numeri- 
cal treatment over the usual steady-state approach 
is that such codes naturally implement the time- 
dependent nature of the SSC process, which is nor- 
mally forgotten in many theoretical investigations 
(e.g. Moderski et al. 2005, Nakar et al. 2009) and 
modeling attempts (e.g. Ghisellini et al. 2009, 
Aleksic et al. 2012). Imagine an electron popu- 
lation that emits synchrotron radiation and then 
scatters this self-made radiation up to 7-rays. The 
electrons, therefore, lose energy, which means that 



the emitted synchrotron emission will also be less 
energetic, as will be the SSC emission. This im- 
plies that the cooling rate will also become weaker 
over time, i.e. the cooling rate is time-dependent. 
Schlickeiser (2009) gave an analytical expression 
for the time-dependent SSC cooling rate, which is 
proportional to an energy integral over the elec- 
tron distribution function itself. 

Schlickeiser, Bottcher & Menzler (2010, here- 
after referred to as SBM) combined the new SSC 
cooling term with the well known synchrotron 
cooling term, in order to give a more realistic 
treatment. As mentioned before, the SSC cool- 
ing becomes weaker over time, which means that 
after some time it will be weaker than the stan- 
dard linear cooling terms, such as synchrotron 
cooling, since they do not depend on time. They 
also calculated the synchrotron SED with the re- 
markable result that the synchrotron component 
exhibits a broken power-law behaviour. Interest- 
ingly, this feature is independent of the electron 
injection function. While SBM used a ^-like in- 
jection, Zacharias & Schlickeiser (2010) performed 
the same analytical calculation using a power-law 
injection. Apart from the high-energy end of the 
synchrotron component the broken power-law be- 
haviour is the same in both cases. This is due 
to the fact that any extended form of the elec- 
tron distribution is quickly quenched into a (5-like 
structure, as long as no reacceleration is taken into 
account. Therefore, the (5-approach is sort of a late 
time limit for any extended injection. 

Following-up on the results described above, 
Zacharias & Schlickeiser (2012, hereafter re- 
ferred to as ZS) calculated analytically for the 
(5-approach the emerging SSC SED. They obtained 
the interesting result that the time-dependent SSC 
cooling leads, in fact, to a dominating inverse 
Compton (IC) component without the need for 
rather extreme parameter settings. Additionally, 
they found that the SSC component also exhibits a 
broken power-law, which may also be independent 
of the injection formj^ 

As we already discussed above, the debate 



Wc note that the type of breaks in the SED that we can 
naturally account for with our model, are usually described 
with rather complicated electron distributions, requiring 
sometimes multiple spectral breaks in the electron source 
energy distribution with practically no theoretical justifi- 
cation (e.g. Abdo et al. 2011a for Mrk 501). 
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whether SSC, EC or both play an important role 
especially in FSRQs is not yet settled. Since we 
showed in ZS that a dominating SSC component 
is easily achievable, it is just straightforward to 
include the EC scenario in our approach. There- 
fore, we will extend the loss rate of SBM with the 
additional contribution of the external Compton 
losses. This will be done in section [21 In section 
[3] we will calculate the resulting intensity and flu- 
ence spectrum of the EC component, where the 
fluence is the time-integrated, i.e. averaged, in- 
tensity. The lengthy details of these calculations 
can be found in appendix |B] and [Cj Transform- 
ing the results into the form of an SED will be 
done in section ^ where we will also summarise 
the results of ZS for the sake of completeness. At 
the end of this section we will give some example 
SEDs with all contributions of synchrotron, SSC 
and EC, and will discuss the results. In section 
[5] we will summarise the results and present our 
conclusions. 

2. The extended loss rate 

Since we intend to calculate the spectrum due 
to external Compton emission, we have to include 
this type of energy transfer between electrons and 
photons in the energy loss rate of the relativis- 
tic electrons. According to Dermer & Schlickeiser 
(1993) the pitch-angle averaged loss rate of elec- 
trons in an external photon field that is isotropi- 
cally distributed in the lab frame, is 
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with the speed of light c = 3 x 10 cm/s, the 
Thomson cross section ctt = 6.65 x 10^^^ cm^, 
the energy of an electron at rest mgC^ — 8.14 x 
10^^ erg, the electron Lorentz factor 7, and the 
Lorentz factor of the radiating plasma blob Ff,. 

We assume that the energy density of the ex- 
ternal photons is isotropic in the galactic (primed) 
frame, and is given by 



AttR'Ic' 



(2) 



where L'^^ is the luminosity of the accretion disk 
surrounding the central supermassive black hole. 



Tsc is the scattering depth of the ambient medium 
scattering the accretion disk photons, and R'^^ is 
the radius up to where the accretion disk photons 
are scattered by the ambient medium. 

Adding the external cooling term to the syn- 
chrotron and SSC cooling terms, we obtain the 
complete electron loss rate: 

\itot\ = \isyn\ + \iec\ + \issc\ 



"7 -f Ao7 / d77 n(7,i) 



9mpC^ 



= Do{l + lec)-f^ + Aoj^ J d7 72^(7, t), 



with Do = AcaTUB/iim^c^) = 1.29 x lO^^fo^g-i^ 
the magnetic energy density ub — S^/Stt, and the 
magnetic field B = 6 Gauss. Schlickeiser (2009) 
gives the constant Aq = So-tCiPoReq/ {nieC^) = 
1.15 X 10^^^i?i56^ cm^^s^^, which was obtained 
during the derivation of the time-dependent SSC 
cooling term. The parameters are given as Pq = 
2 X lO^^erg-is-i, eq = 1-856 x lO^^o^g^.g^ 
Ci = 0.684. We scaled the radius of the spherical 
emission blob as i? = 10^^i?i5 cm. 

The non-linearity of the SSC cooling manifests 
itself in the integral over the volume-averaged elec- 
tron distribution 71(7, i). The factor 
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= 0.09- 
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(4) 



shows the relative strength of external to syn- 
chrotron cooling. For l^c ^ 1 the linear cool- 
ing is dominated by the external photons, while 
for lec ^ 1 the synchrotron process mainly oper- 
ates. Apart from R'^^, which is scaled as R'^^ — 
3.08 X lO^^R'p^ cm, we scaled the quantities in 
equation Q as Q = ^O'^Qx in their respective cgs- 
units, and F^ = lOFfc^i. 

We can now formulate the differential equation 
describing the competition between the injection 
of ultrarelativistic particles with the source func- 
tion S'(7, t) and the energy losses as described by 
equation (|3| inside the spherical emission region 
(Kardashev 1962): 

dn{j,t) d 



dt 



(?7 



[\itotHi,t)]^S{^,t). (5) 



(3) 
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In order to keep the problem simple, we assume a 
monoenergetic instantaneous injection 5(7, t) = 
Qo S il — lo) ^{t)- Inserting equation ^ into 
equation ([5| we obtain 



dn 
~dt 



Do (1 + ?ec) + ^0 



d7 7^71(7, 
= qoH'y- 7o) S (t) 



which apart from the factor (1 + l^c) equals the 
differential equation that was solved by SBM. We 
can, therefore, use their solution. The only differ- 
ence is that we have to insert the factor (1 + lec) 
wherever they have a Dq . The solution is given in 
cases of the injection parameter a, which is pro- 
portional to the ratio of the nonlinear to linear 
cooling at time of injection. We will discuss its 
implications in section [STT] 

For dominating linear cooling (a ^ 1), we get 
the electron distribution 

70 



n(7,t,a < 1) = go (5 7 - 



l + i:'o(l + /ec)7ot 



(7) 

If initially the non-linear cooling dominates 
(a S> 1), we obtain 



X (5 ( 7 — 



n(7,t,a> 1) = qoH[tc-t] 

To 

(l + 3i^o(l + U7oa2t)i/3 

which is valid for times 

- 1 



t < tr 



3Do{l + lechoa^ 



(8) 



(9) 



For late times the linear cooling takes over and the 
distribution function is described by a modified 
linear solution: 



n(7, t,a:^ I) ^ Qq H[t - t J 
70 \ 



(10) 



^ 1^ + Do{l + lechot J 

2.1. The injection parameter a 

We have written the solutions of the differential 
equation ^ dependent on the parameter a, which 
we intend to discuss in greater detail, now. 

We call it the injection parameter, and it is de- 
fined as the square-root of the ratio of the non- 
linear cooling term to the linear cooling term at 



time of injection, i.e. 

l7s.c(t-0)| 



a 



^09070 



|7syn| + |7ec| Doil + lec) 
N501I 



= 2.13 X 10" 



RUl+le 



(11) 



where N denotes the number of radiating parti- 
/|Cles, and we applied the same scaling law for the 
parameters as above. 

For a ^ 1 we see that the linear cooling dom- 
inates, resulting in the solution ([t]). If a 3> 1, 
the non-linear cooling at least initially dominates, 
giving the solutions (|8| and ( [To| ) for early and late 
times, respectively. The time tc marks the transi- 
tion from non-linear to linear cooling. 

Comparing the above given a to the injection 
parameter obtained by SBM (where EC losses 
were neglected), we find that 



CtSBM 



(12) 



(i + U^/"' 

where ass a/ = ^09070/^0- 

This implies that the contribution of the exter- 
nal photons lowers the possibility for non-linear 
cooling. Solving equation (111 for the electron 
density we obtain 



Do{l + lec) 2^ A 

^o7o ^o7o 



2 

a 



2 "^SBM 



QO.SBM- 



(13) 

This also demonstrates the afore mentioned fact 
that with external Compton losses included it is 
harder to cool the electrons non-linearly. Equation 



( 13 1 shows that for the same value of a, i.e. the rel- 



ative strength between linear and non-linear cool- 
ing, one needs a higher electron density in the blob 
compared to the case where the external losses are 
neglected. 

3. External Compton fluence 

The intensity due to inverse Compton collisions 
of electrons with external photons is given by 

Iec{es,t) = Rjec{es,t), (14) 

with the emissivity (see appendix [A| 

CO 00 

u(e) /" , / s / s 
d7 7i(7,t)cr(e5,e,7), 

(15) 
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where e is the normahzed target photon energy 
in units of the electron rest mass rrieC^, is the 
normahzed scattered photon energy, u{e) is the 
target photon density, and 



Sctt 

467^ 



G(g,r) 



(16) 



being the Klein- Nishina cross-section (Blumenthal 
& Gould 1970) with 



G(g,r) = Go{q) 



rV(i-9) 



Goiq) 

r 



2(1 + rg) 

2qlnq + l + q- 2q^ , 
467, 



r(7- 



(17) 

(18) 
(19) 

(20) 



Thus, we obtain for equation (15) 

oo oo 

u{e) 



Idtt 



de- 



d7 



^G(,,r). 

T 

(21) 



Here 



1 



1 



1 



(22) 



denotes the minimum Lorentz factor for the elec- 
trons, below which the electrons would gain energy 
from the photons. 

The fluence is the time integrated intensity 
spectrum, giving an average of the variability in 
all bands, and also incorporating that observation 
times can be much longer than the typical flare 
duration: 



dtl{e„t). 



(23) 



The intention of this paper is to calculate an- 
alytically the complete SED of blazars. We will 
therefore stick to a rather simple approach and 
use the simplest approach possible for the exter- 
nal photon density in the comoving frame: 



"(e) = o^bKc (5 (e - Cec) : 



(24) 



where Cec is the normalized energy of the target 
photons in the comoving frame. Using the electron 
densities Q, ([S]), and (10) we can calculate the 
intensity, and afterwards the fluence. 



3.1. Small injection parameter, a ^ I 



Using equation ([7| in equation ( 14 ) we find for 
the case of a <C 1: 



/ec(es,r) = /oe.(l + r)^GU(r),r(T),e = eec,7 



70 
1 + r 



X H[t] H [70 - 7™„(e = eec)(l + r)] , (25) 

after performing the simple integrations of the 6- 
functions, and substituting r = Do{l + lec)lot- 
Here we defined Iq = RcaTqa^lu'^d {^''^^Ido) ^ 



G g(r),r(T),e = eec,7 = 



7o 



= Go(g) 
r(r) 

'^mini.^ — ^ec) 



l + r 
e^(l-g)(l + r)^ 
270(70 -e«(l + T))^ 



= 4e, 



70 
l + r' 



4eec7o(7o - es(l + t)) 



1 + Jl 



(26) 

(27) 
,(28) 
(29) 



Inserting equation (25 1 into equation (23) we 



obtain after some calculations (see appendix [B]) 
the final expression for the external Compton flu- 
ence in the case of dominating linear cooling: 



F{es) = Fo7o'4/' 



-1/2 



(1 + eece.)3/2 



^cut 



(30) 



where Fa = /o/(Do(l + ^ec)7o), and e^t = 
4eec7o/(l + 4eec7o)- 

In the Thomson-limit Tq = 4eec7o ^ 1 the 
spectrum cuts off at w ro7o ^ 79. The 
Thomson- limit also implies that e^^ ^ (47o)~^- 
Hence, there is no break in the spectrum, and the 



denominator in equation ( 30 ) equals unity. 



In the Klein-Nishina-limit Fq ^ 1 we find Sec ^ 
(47o)~^, and the cut-off at w 79. Thus, the 
spectrum breaks at eg — e^Tc^ . 

3.2. Large injection parameter, a ^ 1 

3.2.1. Early time limit, t <t(. 



For t < tc we use equation ^ in equation (14 1, 
as well as the same substitution for t as above, and 
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obtain after solving the simple integrations 

Iecies,r) = /oe, H [t, - r] (1 + SaV)^/^ 
xG ( g(T),r(T),e = eec,7 = 



70 



(l + 3a2T)i/3^ 
X H[t] H k - 7mm (e = eec)(l + 3aV)i/3l , (31) 



with Tc = {a^ - l)/3a^, 
G I g(r),r(T),e = eec,7 



70 



Go(g) 



(l + 3a2r)i/3 
e2(l-g)(l + 3aV)2/3 
27o(7o-e.(l + 3a2r)i/3) 



, (32) 



r(r) = 4e, 
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(l + 3a2r)i/3' 

e,(l + 3aV)2/3 
4eec7o(7o-es(l + 3«2T)i/3) 



(33) 
(34) 



Using equation \2'M with equation (31 1 the flu 



ence for the early time limit can be calculated. 
The lengthy details can be found in appendix [C| 
The solution depends on the external photon en- 
ergy Eec, giving three different cases. 
For e^^i < 2(%/2 - 1)7^ we have 



2(y2-l)7B 

X ( 1 - ^ ) , (35) 

^CUt , 



while for 2(\/2 - 1)75 < < 2(\/2- 1)70 we get 



F{es) = 



5/2 



(1 + eecEs) 



5/2 



X 1- 



^cut 



(36) 



The last part is for 2(-\/2 — 1)70 < e^^ and becomes 



3.2.2. Late time limit, t >tc 

As before, we can find the intensity for t > t^, 



if we insert equation (10) into equation (14 1. As 



in the previous cases we exchange t, and get 
Iecies,r) = Joe, H [t - rj (^^^ + rf 



3a2 



xG |^q(T),r(r),e = eec,7 = 

r 1 + 

X iJ 



70 



3q2 



, ,,l + 2a''' 
70 - 7m*"(e = eec)( - H t) 



3a2 



(38) 



with 



G (7(T),r(r),e = tec, 7 



7o 
1 + r 



Go(<z) 



g^(i-g)(45^ + r)^ 

27o(7o-e«(ii|?^+r))' 



9(t) 



70 



3a2 



e.(l + r)2 



4eec7o(7o - es(l + t)) ' 



(39) 

(40) 
(41) 



Inserting equation ( 38 1 into equation ( 23 ) one 



can obtain the fluence in the late time limit. For 
the details we refer the reader to appendix [C] 
again. The results depend also on the external 
photon energy, giving two cases this time. 
For < 2(^/2 - 1)75 we find 



IB I 



FoYo^ 



3 3/2_ 



^1/2 



(l + eeces)3/2 



1 - 



(42) 



while for e-,^ > 2{V2 - 1)75 we obtain a single 
power-law in the form 



(43) 



F{es) = 



5/2 



2(\/2-1)7b, 

x(l-^ 



(37) 



Here we used jb = 7o/ct, and e^^ = 4eec7B/(l 

4eec7B)- 



3.2.3. Total fluence in the case a ^ 1 

Combining the results of the previous sections 
we obtain the total fluence for the inverse Comp- 
ton component due to interactions with the ambi- 
ent radiation field in the case that the electrons are 
at first cooled non-linearly. We have three differ- 
ent cases depending on the value of the normalized 
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external photon energy e^c- 
For < 2(%/2 - 1)7B we find 



-1/2 



xfl-— ). (44) 

V Cent / 

Secondly, if 2(^2- 1)75 < < 2(^/2 - 1)70 the 
fluence becomes 



-1/2 



1+ ^) (l+eecEs) 

e 



5/2 



X ( 1 - — ) . (45) 

^CUt , 



Lastly, we obtain for 2("s/2 — 1)70 < e^^ 



-1/2 



5/2 



2(V2-1)70 

X ( 1 - ^ ) . (46) 

^cut , 



These are all broken power-laws, where at least 
one break (e-y^) depends strongly on the value of 
a, and therefore on the non-linear cooling. 

We note that the first case corresponds to the 
extreme Klein-Nishina- limit (Fq 1), the second 
one to the mild Klein-Nishina-limit (Fq > 1), and 
the last one to the Thomson-limit (Fq !)• 

4. The complete SED 

Using the results of the last section and of 
we are now in a position to present the complete 
spectral energy distribution in a combined picture 
of SSC end EC radiation. The inverse Compton 
component will, therefore, be the sum of the SSC 
and EC contributions. If these contributions are 
comparable on some scales the resulting spectrum 
will deviate from pure power-laws. 

In order to give the results in a manner that can 
be easily compared to data, we will give the SEDs 
depending on the frequency v' = ^^^5e, which is 
already transformed to the frame of rest of the 



^Some of the results of ZS have printing errors, which we 
correct here. 



host galaxy with the Doppler factor 8 = (F;,(l — 
PrbfJ'obs))~^ , where /io6s is the cosine of the angle 
between the jet and the line of sight, and = 
Vb/c is the normalized speed of the plasma blob. 
The SED is then given by the fluence multiplied 
with the frequency /' — ncRv' F' (v') in units of 
erg/s with the transformed fluence F' — 5^F. 

We should note that we have to adapt the re- 
sults of ZS to the case discussed here. That is, we 
have to include the addition of the external Comp- 
ton cooling to the linear term. This can be done by 
replacing any "Dq" of ZS with " Dq{1 + lee)". In- 
terestingly, the synchrotron component will not be 
affected by this replacement, while the SSC com- 
ponent gains a factor (1 + lee)- 

Below we will first present the theoretical SEDs 
and afterwards give a brief discussion of the re- 
sults. 

4.1. Synchrotron SED 

From equation (ZS-69) we obtain the syn- 
chrotron SED for the case a <^ 1, which is 



a / V 



74 \ l^syn 



1/2 



erg/s, 
(47) 



with Vsyn = 4.1 • 10"567| Hz. 

The maximum value of the synchrotron SED, 

2 

= 2.4 • 1038i?i554^ erg/s, (48) 



74 



is attained at 



<.ma. = ^^syn = 2.1 • 10^^5674 Hz. (49) 



Using equation (ZS-82) we can write the syn- 
chrotron SED for a ^ 1 as 



= 5.6 • W^Ri^S^ 



1/2 



74 1 + 



It peaks at 



erg/s. 
(50) 



ve = 2.9 • 1014(5^ Hz (51) 



with the maximum value 



a 



= 2.4- 103«i?i554- erg/s. (52) 

74 
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4.2. Synchrotron self-Compton SED 

For a ^ 1 ZS found two different versions of the 
SSC SED depending on the Klein-Nishina parame- 
ter K = O.I3667I. In the Thomson Hmit {K «; 1) 
we use equation (ZS-72) and get 



fssciK < 1) = 2.8 • 10-^''i?i5<5^(l + lec) 



a 
74 



VT ) 



e-'^V'^T erg/s, (53) 



while for the Klein-Nishina limit {K ^ 1) with 
equation (ZS-73) we have 



fssciK » 1) = 2.1 • 104°i?i5^^(l + lec) 
1 



ML 



7/4 



7/12' 



erg/s, (54) 



with the constants vt = 1-6 • l(fi^8h^\ Hz, the 
break frequency vb = 2.Z-W^5h-^/^ Hz (Schlick- 
eiser & Roken 2008), and v^^ = 1.2 • 10^4^74 Hz. 
The maximum frequencies 



^SSC,ma.iK « 1) = ^l^T = 1-2 • lO^^^^T^ Hz, 



(55) 



and 



'^SSCmaAK » 1) = ^tys = 1.7 • 102*56-1/3 Hz 



imply the maximum values 



(56) 



fssc,ma.{K « 1) = 1.0-10^^ RuS\l+lec)^ erg/s, 

(57) 



74 



and 



fsscma^K » 1) = 6.3-10^^ Ri^S^l+hc)"^ erg/s 



hi 



respectively. 



(58) 



For the SSC SEDs in the case a » 1 we have 
to discuss three different cases depending on the 
Klein-Nishina parameter K, and we begin with the 
Thomson limit {K <C 1) using equation (ZS-85): 



fssciK « 1) = 2.8 • 1039i?i5(54(l + /,,) 



a 



74 



, \3/4 



The maximum value 



4 



erg/s. (59) 



fsSC,ma.iK « 1) = 1.6-10'^ Rl5S\l+lec)- erg/s 

(60) 



74 



is attained at 



'''sscmu.iK « 1) = -^i^T = 4.0 • lO^^SlTft Hz. 

(61) 

For the mild Klein-Nishina limit (1 <C if <C a^) 
equation (ZS-88) yields 



fssci^ < if < a^) = 1.5 • 1039i?i554(i + 



74 



^x3/4 



(l+^) (l + ^) 
This triple power-law peaks at 



13/3' 



erg/s. (62) 



'^'sSCma.i^ « K « a^) « ^l^B = LQ-lO^^^fc- V3 Hz, 

(63) 

and reaches a maximum value of 



f'sSCmaA^ « K « a') « 1.2.1039i?i5^'(l+«ec),T7^ Crg/s 



a 

'6V37I 
(64) 

The last case is the extreme Klein-Nishina limit 
(K < if). We obtain the SED from equation 
(ZS-91), resulting in 



fssci^ < a' « if) = 1.1 • 104°i?i5r(l + kc) 



a 



, \3/4 



X 



n-F ;^^^) (1 I 



1- 



70 



13/3' 



erg/s. (65) 



The SED peaks with a maximum value 
fsscmaA^ « « ^) ~ 3.4.1039i?i5<5'(l+«ec)^ erg/s 



&7l 
(66) 



at a peak frequency 



'^'sscmaA^ « « ^) « 4^^^ = 1.7.10^*56-^3 Hz. 



(67) 
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4.3. External Compton SED 

In the case of a <C 1 the SED can be calculated 



from equation ( 30 1 , becoming 



f,,iiy') = 4.0 X m'^Ri^lec^- 
— , 3/, 1 erg/s, 

^ , Vent) 



(68) 



where i^br ~ 1-23 • Vd^^bt^}^ Hz, and v^ut — 
1.23 • W^btcut Hz. In the Thomson-Umit (cec < 
(470 )~^) the maximum frequency 



ec,maa; 



1.6 X 102°(5eec7o Hz 
leads to a peak value of 

(-4 9 

a 

74 



(69) 



f 

J ec 



4.4 X 103«i?i5?ec^^ erg/s. (70) 



In the Klein-Nishina- limit {eec ^ (47o) ^) the 
maximum value 



ec.max ~ 1-6 X W^'^Ri^lec^ 

74 ^ec 



is attained at 



erg/s (71) 



(72) 



Now, we convert the results of equations ( 44 1 , 



(45), and (46) into SEDs, which are the cases for 



a 3> 1, yielding for the extreme Klein-Nishina- 

limit Vbr < V^gr 

/^,(z/')=4.0-1034i?i5Z,,4^ 

(^) y 



1 - 



Vcut 



with 



1.02 • 10^4,5^ Hz, and v-y 



erg/s, (73) 



1.23 



lO'^'^Se^g Hz. It peaks at 



1 



Vbr = 6.2 • 10^^ 5e-^ Hz 



(74) 



with a maximum value of 



15''ec o 
74 



erg/s. (75) 



In the mild Klein-Nishina-limit v^^r < v^r < 
ly-fgr, with v^gr = 1-02 • 10^^(574 Hz, we obtain 



f'Jv')^ 4.0- 10-'^ Rl^lec^- 

74^60 

(^) ^ 



s 5/2 I 1 - — ) e^'g/S- (76) 



The maximum value 



/ecna. « 2.0 • 10^^i?15^ec ^37^ erg/s. (77) 

74 

is attained at 



™„„ w = 1.23 • lO^'^S— Hz (78) 



ec,ma2: 7-B 



Finally, we have the Thomson- limit Vy^r < '^far, 
where the SED becomes 



- 2 

74 fee 



/^,(i.') = 4.0 • lO^^i^is^, 



, \l/2 



1+7^1 1 



1 - 



erg/s. (79) 



It peaks at 



ec,max 7-B 



4.92 . 10^85^ Hz (80) 



(note that we use a different approximation for 
i^jg, that is for e^g, compared to equation (78)) 
reaching a maximum value of 



J e 



ec.max 



r4 

1.9- lO^s^is^ec— erg/s. (81) 

74 



4.4. Discussion 

In figure [l] we present several model SED for 
different sets of parameters. Since the number 
of free parameters is fairly large a strict param- 
eter study is beyond the scope of this paper, al- 
though we will discuss some aspects below. There- 
fore, we use some global parameters, which are 
unchanged through all figures, namely S = 10, 
and Th^i = L'^q = r_2 = Rpc = ^15 = 1- The 
normalized external photon energy is also fixed at 



9 




Fig. 1. — Model SEDs for several cases of a, lec, and 79. Specific parameters are given at the top of each 
plot, while the global parameters are given in the text. Top row: a <C 1 in the Thomson limit for both cases 
of lec- Middle row: a 3> 1 in the Thomson limit for both cases of Zee- Bottom row: a 3> 1 and lec ^ 1 in 
the mild and extreme Klein-Nishina regime, respectively. 
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Cec = 10~^. This leaves us with three free pa- 
rameters: the injection parameter a, the relative 
strength of the external Compton cooling l^^y £^nd 
the initial electron Lorentz factor 70. These pa- 
rameters are given at the top of the plots for each 
case. 

Since we use lec as a free parameter and fixed 
almost all parameters it contains, we get the mag- 
netic field strength 6 as a dependent variable. In- 
creasing lec results in a decreasing magnetic field. 
This choice also demonstrates that the figures 
given as examples here only cover a very narrow 
range of possible realizations. 

We note that for the top panels of figure [T] we 
use equation (68) for the EC component, while we 



use equation (79) in the middle plots and equa- 



tions (76), and (73) for the lower plots, respec- 



tively. This is due to our choice of parameters. 

The first obvious result is that for a <C 1 the 
inverse Compton component is significantly influ- 
enced, if not dominated by the EC component. 
For a ^ 1 the main contribution comes from the 
SSC component, at least in the Thomson- and 
mild Klein-Nishina-limit. This is expected, since 
a controls the main cooling mechanism, which 
should also influence the SED. 

Similarly, this is also correct if one considers the 
linear coolings only, where the relative strength is 
controlled by l^c- For example, in the top plots 
of figure [1] where we use equation (68), one can 
see that for l^c ^ 1 the synchrotron component 
dominates the EC component, while for l^c ^ 1 it 
is the other way around. Since lec = |7ec|/|7si/n|, 
one could expect that the ratio of the peak val- 
ues Rec,s = f'ec.max/f's.rnax « 'ec- This is Certainly 



true for the Thomson limit of equation ( 68 ) , as is 



shown in equation (D2) of appendix D However, 



due to our choice of parameters the top panels 
of figure [T] use the Klein-Nishina-limit of equation 
(68), where equation (D3) demonstrates that the 



ratio is more involved and depends also on 70 , and 
Cec- Using the parameters of the plots in the top 
row of figure [l] we find that the ratios become 
Rcc,s ~ 0.07 and Rec.s ~ 66, respectively, consis- 
tent with the plots. 

In appendix |D] we also show the ratios be- 
tween the maximum values of the SSC and the 
EC components. As one can see, every Rssc,ec oc 
{1 + lec) /he with the consequence that the ratio is 



independent of lec for lec ^ 1- Although the de- 
pendence on other parameters is rather involved 
and changes from case to case, we can say that a 
lower value of the external photon energy Cec fa- 
vors the EC component, while a higher value of the 
injection parameter a favors the SSC component, 
as one could expect. 

On the other hand, if lec ^ Ij we find Rssc.ec c< 
lec ■ This implies that the EC component is much 
reduced compared to the SSC component. This 
is also expected, since in this case the EC cooling 
is the weakest operating mechanism, if the other 
parameters are chosen accordingly. It is, of course, 
still possible that SSC could be lower than EC, 
since the ratios in appendix [D] can always become 
smaller than unity. 

As one can see in the lower panels of figure [l] 
we achieve similar results in the mild and extreme 
Klein-Nishina-regime, respectively. Especially in 
the latter case we have an example that the in- 
clusion of the EC cooling has significant effects on 
the emerging spectrum, giving a different power- 
law and peak value compared to the case where 
either EC or SSC is neglected. 

As stated above, in figure [T] we just show some 
model plots with only a very limited range of appli- 
cation, since we fixed most of the free parameters. 
Some of them, however, might have a significant 
impact on the resulting SED, like the radius R of 
the plasma blob, which we set to i?i5 = 1. This 
value corresponds to a light crossing time scale of 
tie = i?/((5c) « 0.9h (with 5 = 10), which is usu- 
ally equaled with the variability time scale t^ar of 
the source. Naturally, a broad range of variabil- 
ity time scales is observed in blazars, ranging from 
minutes (as in the cases of PKS 2155-304, Aharo- 
nian et al. 2007, and PKS 1222-^216, Tavecchio 
et al. 2011) to months (e.g. Abdo et al. 2010), 
indicating a large range of possible radii for the 
emission blobs. 

Increasing the radius of our plasma blob would 
result in a decreasing value of the injection param- 
eter, since a oc R~^ according to equation (11). 



This implies that the SSC component would be- 
come less important, leading to a dominating EC 
fiux in the high energy regime. 

On the other hand, long variability time scales 
can also be the result of inefficient cooling or other 
factors with long time scales, such as probably a 
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rotating jet (e.g. Agudo et al. 2012). Our sim- 
ple model assumes an instantaneous and effective 
cooling of the electrons far away from steady-state, 
favoring short variability time scales, at least if one 
is interested in the effects of the non-linear SSC 
cooling. Schlickciser (2009), as well as Zacharias 
& Schlickeiser (2010) showed that the non-linear 
process begins orders of magnitudes earlier to cool 
the electrons than the standard linear synchrotron 
cooling, implying shorter variability time scales. 
However, a true analyses of the variability should 
be performed with the help of light curves, which 
we aim to present in another paper. 

We should also note that our discussion adopts 
a point-like emission zone, which means that we 
neglect additional effects from photon retardation 
variability processes from the finite size of the 
emission zone. For an extended analysis of these 
effects we refer to the discussion in Eichmann et 
al. (2010, 2012). 

Another parameter, which might easily have 
different values, is the Doppler factor 6. In the 
literature it usually covers a range 10 — 30, but 
extreme values of 50 and higher have been used 
to model blazar spectra (e.g. Begelman, Fabian 
& Rees 2008), although the latter are mostly re- 
garded as unrealistic. 

For our plots we used S = 10. Increasing that 
value would not change the overall appearance of 
the SED, since in our simple model the dependen- 
cies on S are the same in all cases, however the 
peak luminosities would increase significantly. In 
the plots in figure [l] we achieve 7- ray luminosities 
roughly between 10'^^ erg/s and 10^^ erg/s. The 
former value is rather low for blazars. However, 
the values of the peak luminosity depend strongly 
on the parameters, which might differ from the 
ones used to create the plots. 

As an example: With S = 50, a = 50, and 
lec = 100 one would arrive at a SSC luminosity 
of roughly 10^" erg/s. This is close to the high- 
est recorded luminosity of blazars, achieved by 3C 
454.3 with a peak luminosity of 2.1 x 10^" erg/s 
(Abdo, et al. 2011b). In principle, the above 
stated values for a and l^c could be achieved by 
setting A^5o « 120, and b = 0.14. The high number 
of electrons in the source may not be that unre- 
alistic, since a large outburst requires something 
extreme going on, such as a large mass load into 



Table 1: Energetic estimates of the plots in figure 

m 



Plot 


lie [ erg cm ^] 


Pe [ erg/s] 


Ue/UB 


TL 


9.1 X 10-* 


9.2 X 10^'' 


3.0 X 10-^ 


TR 


9.1 X 10^2 


8.5 X 10''! 


2.4 X 10^ 


ML 


9.1 


8.9 X 10''-'' 


2.7 X 10^ 


MR 


9.1 X 10^ 


8.5 X 10*5 


2.5 X 10^ 


LL 


91 


8.5 X 10*4 


2.5 X 10^ 


LR 


9.1 


8.5 X lO'*^ 


2.5 X 10^ 



Note.— TL: top left, TR: top right, ML: middle left, 
MR: middle right, LL: lower left, LR: lower right 



quency of the SSC component would be at roughly 
3 X 10^'^ Hz, which is reasonable (Abdo, et al. 
2011b). Using equation (51) we obtain a syn- 
chrotron peak frequency roughly at 8 x 10^"'^ Hz, 



which is in the far infrared regime. Still, a Doppler 
factor of 50 is quite high, severely questioning the 
validity of the results of our simple model. Only a 
detailed spectral analysis, potentially also includ- 
ing a discussion of variability aspects (see above), 
can decide about the validity. However, such an 
analysis is beyond the scope of this paper. 

Another important aspect is the energy trans- 
ported and dissipated by the content of the jet. 
The initial comoving energy density of the elec- 



trons is given by Ue 
( 13 1, we obtain 



qo'-foTTieC . Using equation 



(l + lecW -3 

We = 0.09 erg cm . 

-R1574 



(82) 



The power transported by these electrons is given 
by Pe = TiR^clST.Tlue, yielding with Pr, ~ 1 

Pe = 8.4 X 10«r2 ^rJ1±^^ erg/s. (83) 

74 

A widely used measure to quantify if the jet is 
dominated by particles or the magnetic field, is 
the equipartition parameter 



(84) 



which can be rewritten with the help of equation 
Q, giving 



the jet. According to equation (61) the peak fre- 



24.4- 



(1 + lec)leca'^ 



UB 



(85) 
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Table [T] lists the values of these constraints for 
the plots given in figure [T] Depending on the used 
parameters we obtain a wide range of possible re- 
alizations. 

A first result is that we are far away from 
equipartition with all our set-ups, and only the 
case with the lowest energies involved shows a 
dominating magnetic field. This is reasonable, 
since the number of electrons in this case is also 
at a minimum in our models. This is also the only 
case, where the synchrotron peak dominates the 
IC peak. However, we should note that this might 
be an artifact of our parameter settings, i.e. the 
inverse dependence of b on Zeci while we left all 
other parameters in l^c unchanged. Keeping b at 
some value and changing some other parameters 
in Zee might give a different picture. 

The energy density and the transferred 
power Pe of the electron naturally increase with 
the number of particles in the blob. Comparing 
our values of Pe with, e.g., Ghisellini et al. (2009), 
we find that our plots cover a range that is usu- 
ally below obtained values for FSRQs, which is 
again due to our choice of parameters. The case 
with a — 10, lec — 100, and 74 = 1 is, however, 
exactly in the range obtained by these authors. 
Higher values, as we used for the discussion about 
3C 454.3 above, are also possible. So, from the 
energetics point of view our model is definitely 
supported. 

5. Summary and conclusions 

In this paper we aimed to extend the earlier 
work of Zacharias & Schlickeiser (2012, ZS) re- 
garding the importance of time-dependent syn- 
chrotron self-Compton cooling scenarios for mod- 
eling the spectral energy distributions of blazars. 

In ZS we were able to give an analytical SED for 
a pure synchrotron/SSC picture. However, espe- 
cially in the context of fiat spectrum radio quasars, 
a preference is given in the literature to external 
Compton scenarios, where the high energy compo- 
nent is modeled using external photon sources that 
are inverse Compton scattered by the jet electrons. 
We considered this effect here by introducing the 
external Compton cooling term into the differ- 
ential equation describing the time- and energy- 
dependent behaviour of the electron distribution 
function. We used a rather simple model, where 



we assumed the external photons to be isotropi- 
cally distributed in the frame of the blob. This is 
justified by our intention to show in an analytical 
demonstration the possibilities and differences of 
our approach compared to the usual steady-state 
approach. More realistic scenarios can mostly be 
treated only numerically, since they contain more 
sophisticated assumptions and effects. 

We then calculated the resulting external 
Compton intensity, the fluence and eventually the 
SED. Adapting the results of ZS regarding the syn- 
chrotron and the SSC SED, we have now obtained 
a complete analytical prediction of the appearance 
of blazar SEDs, if the electrons are cooled by the 
non-linear, time-dependent SSC channel in com- 
bination with the linear synchrotron and external 
Compton channels. 

The main addition to the conclusions of ZS is 
that it critically depends on the parameters of the 
source, which type of cooling dominates. For only 
weak EC cooling the results of ZS are recovered, 
while in the other case the shape of the high- 
energy component can change significantly. In ap- 
pendix[D]we calculated the ratios of the maximum 
values of each component with the other compo- 
nents, since they mark the dominance of one com- 
ponent over the others. These ratios especially 
show the importance of the choice of the specific 
values of the parameters, and demonstrate that a 
broad range of realizations is possible. 

We also note that the nonlinear cooling has sig- 
nificant effects on the shape of the spectrum of all 
components leading to additional and also unique 
breaks, which strongly depend on a. This high- 
lights once more the severe effects of the time- 
dependent treatment. 

To conclude, we can say that our investiga- 
tion so far has clearly shown the importance of 
the inclusion of the time-dependent nature of 
the SSC cooling term, especially in blazar flar- 
ing states, which are far from steady-state. The 
resulting spectra differ significantly from the usual 
approach, where these effects are not taken into 
account. Secondly, due to the broad range of free 
parameters in the EC model, the SSC component 
should not be dismissed beforehand while model- 
ing blazars. As we were able to show here, the 
SSC component might be as important as the EC 
component, if the correct time-dependent treat- 
ment is applied, which is especially necessary in 
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flaring states. 

In order to quantify the difi^erences between the 
time-dependent and steady-state approach a bit 
further we intend to calculate the resulting light- 
curves and discuss the effects of optical thickness 
in a future work. The former might show in a 
clearer form that the SSC cooling actually acts a 
lot quicker than the usual cooling scenarios, while 
the latter might be important to constrain the free 
parameters, such as the Doppler factor. 

We thank the anonymous referee for the con- 
structive and detailed comments, which helped 
significantly improving the manuscript. 
We acknowledge support from the German Min- 
istry for Education and Research (BMBF) through 
Verbundforschung Astrotcilchenphysik grant 05A11PC1 
and the Deutsche Forschungsgemeinschaft through 
grant Schl 201/23-1. 
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A. The external Compton emissivity 

According to Dermer & Schlickeiser (1993) the general emissivity due to external Compton scatterings is 
given by 



ce 



oo oo 

J de j) dn J dj ^ dl7e(l-/3rcos*) 



(Al) 



where f2, and He mark the direction of the incoming photon and the electron, respectively, and cos 5* 
lifj,e + (1 — fj,'^y^^{l — /^e)^^^ cos ($ — <I>e) is the angle between the incident photon and the electron. 
Assuming complete isotropy we can write u{e,Q) — M(e)/47r, 71(7, fJe) = ^{'j) / An , and (T(es, e, 7, f2, rjg) 



(T(es, e, 7)/47r. Plugging this into equation (All the angle integrations can be easily performed, yielding 



47r 



00 00 



(A2) 







which agrees with equation (15 1 



B. Intensity and fluence for a <C 1 

Using equation ^ in equation (14 1 we find for the case of a ^ 1: 



d7^G(9,r)<5(7- 



70 



1 + Doil + lechot 



(Bl) 



which leads directly to equation (25) with the substitution r — Do(l + lec)lot and the definition Iq 
Rc(JTqQTlu'J{Anel^-ll). 



Integrating equation ( 25 1 with respect to time we yield 

10 I 
Imin 

F{es) = Foe, J dr {1 + t)^G (^q{T),T{T), 



e = eec,7 



70 
1 + r 



where we set Fq = /o/(L'o(l + ^ec)7o)- 
Inverting (/(r) gives 



T{q) ^ 2eec7o9 U/ 1 + 



-1-1. 



This can be substituted into equation (B2l yielding 



F{es) - 4FoeL7oes J 



dq 



1 



1 



^s^ecq 



-1 Giq). 



With the above substitution the function G simphfies to 

G(g) -Go(g) + 263e,,g(l-g). 



(B2) 



(B3) 



(B4) 



(B5) 
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The final substitution x ;= (esCec'?) ^ gives 

(yrr^- 1) 



7^ 



dx 



Gq 



1 



(B6) 



This integral cannot be solved in a closed form. However, we can achieve approximative results giving 
reasonable solutions for most of the parameter space. 

For < one can see that both integration limits are much larger than unity, and so is x. Thus, we 
can approximate the fraction 

(^/^T^-l)' 



-5/2 



X-^y/l + X 



(B7) 



The bracket in equation (B6) is to leading order in this approximation G{x) ~ 1. Now, this can be easily 
integrated, yielding 



72 



^0-73^3/2 -1/2 



1 



3/2' 



(B8) 



For €s > and > 2(-\/2 — 1)70 both integration limits are smaller than unity the fraction in the integral 

,4 



of equation ( B6 ) becomes 



(Vl + a;-l)_^2-4. 



x'^y/lTx 

The leading order of G{x) is here G{x) sa 2/ [esf-ecX^)- The integral, thus, results in 

1 



(B9) 



F(e, >2(V2-l)7o)«4Fo7ile; 



dx 



3,-2 



4eec7o(7o - Es) 



(BIO) 



For intermediate frequencies t^^ < < 2(%/2 — 1)70 we can split the integral into two regimes and use the 
above approximations: 



F{e-J < < 2(V2 - 1)70) « 4Fo7o'e 



3^-2 



dx 



236 



da;x-^/2 



19 

1) 



3^-2 



(Bll) 



Combining the above results we obtain the final expression for the external Compton fluence in the case 
of dominating linear cooling: 



-1/2 



(1 + eeces)^/^ V" ^cut) 



1 - 



(B12) 
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where ecut = 4eec7o/(l + 4eec7o) is the solution of the equation es/(4eec7o(7o ~ ^s)) — 1, yielding equation 
dSOl). 



C. Intensity and fluence for a 3> 1 
C.l. Early time limit, t < tc 

For t < tc we use equation ^ in equation (14) and obtain 



IQtt 



90 



7o 



7 V {l + 3Do{l + lechoaHy/-^ 



US -(CI) 



Using the same definitions of r and Iq one can easily arrive at equation (31 ). Then, the fluence integral can 
be written as 



oo 

F(e,) = Foe, J dr (1 + 3a\)^/^G (^^(t), r(T), e = 



70 



(l + 3a2r)i/3 



H Wr. -t] H 



70 



'^min (^ec) 



(l + 3aV)i/3 



.(C2) 



Examining the Heaviside functions we find that is the upper limit as long as < , while for > e^g 

3 1 

/3q;2. As long as they arc not really needed we will refer to both of 



the upper limit becomes 



1 



them as the respective variable with the subscript up. 
The break energy is given by 



IB 



1 



(C3) 



47b Ee 



where 7b = 7o/a- 

Substituting ^ = (1 + Sa^r)^/^ we obtain 



(C4) 



where fiupi^s < £7^) = and fiupi^s > £7^) = 7o/7n 
Inverting q{fi) yields 



fi = 2eec7o9 \/ 1 



- 1 



(C5) 



Inserting this into equation (C4) with Qupi^s < £73) — £sa^/(4eec7o(7o ~ ^sCt)), and qupi^s > £7^) — 1, we 
find 



^ v^sj - £s 



gup 



dq 



1 



l + — 1) [Go((7) + 2e,eecg(l-g)]. (C6) 



4eec'70{T0-^s) 

Using again the substitution x = {tstecQ)~^ we eventually obtain 



16Fo7o' -4 



a 



2 



da; 



(Tm-i)^ 



Go 



(C7) 
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with the lower hmits Xdowni^s < £73) = 470(70 - aes)/e^Q!^, and Xdownies > £jb) = (esCec) 

Since this integral cannot be solved in a closed form, as well, we will now consider similar approximations 

as in the linear case. Additionally, we will now consider each case of the lower limit in turn. 



Be ginning with, < es < ecut we use the lower hmit Xdowni^s > ^-yg) = (esCec) ^- In case of Cec > 
(2(^2-1)75)-' we find for 2(^2-1)70 < £s < Scut that both limits are lower than unity and we approximate 
the fraction in equation ( |C7[ ) with 

r6. 



x'^Vl + X 



(C8) 



and G{x) becomes to leading order G{x) sa 2/{eseecx'^). Hence, 



F(2(V2-l)7o<es<ec„t) 



I6F075 



dx 



1 



2a 



2 s 



^cut 



(C9) 



For Cyg < Cs < 2(%/2 — 1)70 we see that the upper limit is larger than unity, while the lower limit is still 
smaller than unity. We will, therefore, split the integral. For the part below unity we use the above mentioned 
approximation, while for the part larger than unity we approximate the integrand with x~'' 1"^ . Thus, 



F(e^^ <e, <2(n/2-1)7o) 



9 S 



dx 



(e,eec)-i 



da;x-^/2 



2a2 



(CIO) 



Combining the last two equation we obtain 



Foil ,-4 



2a 



2 s 



1 - 



(Cll) 



For < (2(-\/2 — l)7s) ^ we can make use of equations (C9| and (CIO), while a; ^ 1 and a; « 1, 
respectively, yielding 



F(2(V2- 1)70 <es< e^^t) « ^e^" 



1 - 



f-cuL 



(C12) 



and 



F(e-i <e, <2(V2-l)7o) 



J^o7o'.-4 



2a 



2 s 



(C13) 



If both limits are larger than unity we approximate the integrand with x and obtain 



■*T0(T0-'a) 

— 



I6F075 



dxx-'^/^ 



32-Fo7ogec^ _3/2 
5a2 



1 - 



^cut 



5/2' 



(C14) 
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Combining the equations (C12) - (C14| we can write the final result as 

-3/2 



F{ejg < es < tcut) 



2a2 (i + e^^£^)5/2 



1 



- 



^CUt J 



(C15) 



Now, we turn our attention to the case where < e^^ , which means that the lower limit in the integral 
of equation (C7) becomes Xdown — 47o(7o — ce^s)/('loi^- 

The first thing, we notice, is that the upper limit is always larger than unity. This is also true for the lower 
limit, except for e'^^ < 2(\/2 — 1)75. In the latter case, we obtain with the same approximations as above: 



16i^o7o 



dx 



42Q{jQ-aes)_ 



'I'l 



(C16) 



In the case, where both limits are larger than unity, we find 

b [ts <<^-,b) 5 — 







7^ 



(C17) 



Finally, we can give the fluence in the early time limit covering all energies. However, we have to distinguish 
three cases for the external photon energy. 
For ej^i < 2(\/2 - 1)73 we have 



F{es) - 



while for 2(^2 - 1)75 < e^^^ < 2(%/2 - 1)70 we get 



2{^/2-l)tB 



1 - 



^cut 



(C18) 



Foa^ 



The last part is for 2(\/2 — 1)70 < e'^^ and becomes 



^cut 



(C19) 



He.) 



2(^/2-l)^ 



5/2 



^cut 



(C20) 



Hence, we obtained equations (35) - (37). 
C.2. Late time limit, t > tc 



As before, we can find the intensity for t > t^, if we insert equation (10) into equation (14): 



Idtt 



90 



7o 



^ + Do{l + lechot 



l(.C21) 
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Equation (38) is recovered if one uses the definitions of r and Iq, again. Tlie fluence then follows to be 



10 l + 2o 



F(e,) = Foes 
Substituting y = -^"^^"^ + r we obtain 



dr 



3a2 



7o 



3a^ 



(C22) 



F{es) = FqEs J dyy^G ^q,r,e = eec,7 = 



(C23) 



Inverting q{y) gives y = 2eec7o'7 ( + ^ l) ' which leads to 



F(e,) = 8i^oe^c7oe. 



dg 



1 



[Go{q) + 2ese,,q{l-q)]. (C24) 



Finally, we use again x — (eeces?)^^ and, thus, the fluence becomes 



F(e,) = 8^^o7oe;' / da: 



Go 



1 



1 



(C25) 



As before, we will consider approximative results. For e^,} < 2{\/2 — 1)75 < Es < e-y^ both limits are 
smaller than unity, and we approximate the integrand with {2^eseecx'^)^^ giving 



F{2{V2 - 1)75 < < e^J « SFoj^e-^ [ dx 



(e.eec)-i 



(C26) 



For < e^^ < 2(\/2 — 1)75 < £7^ both limits are larger than unity and the integrand can be written 
a;~^/^, yielding 

e2a2 

F{es < < 2(72 - 1)75) « SFoj^eJ^ J dxx-'^' 



as 



16 

y 



3/2' 



(C27) 



In the intermediate part e^^ < < 2(\/2 — 1)75 < e-y^ we can use both limits to obtain 



F(e-i < < 2{V2 ~ 1)75) « 8i^o7oe 



3,-2 



47o(70-Qes) 
7^ 



dx 



1 



i^o7oeJ • (C28) 
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Collecting terms we find the fluence in the late time limit for e^^ < 2(^/2 — 1)7b 

-1/2 



3.3/2 _ 



1 



For e^^ > 2(%/2 — 1)7b we obtain a single power-law in the form 



(C29) 



(C30) 



Both equations equal equations (42) and (43), respectively. 



D. The ratios of the maximum values 

The ratios of the maximum or peak values are defined as 

f 

•f I. max 



f 



(Dl) 



D.l. External to synchrotron 



In the case a <C 1 we use equations (48) and (70) for the Thomson- limit of the EC component, and obtain 

(D2) 



Rec.s — 1-8 le 



In the Klein-Nishina- limit of the EC component we use equation (71 1, and get 



Rec,s = 6.6 X 10 



-5 'ec 



74 



(D3) 



For a 3> 1 we have to distinguish between the three cases for v^r ■ For the extreme Klein-Nishina-limit we 

(D4) 



use equations (52) and (75), giving 



R^^ ^ = 6.6 X lO^-'le 



74 ee 



For the mild Klein-Nishina-limit we use equation (77), yielding 



R,r , = 8.3 X 10~-^L 



a 



74^6 



1/2 



while we use equation (81) for the Thomson-limit, resulting in 

Rec.s — 0.8 lec- 



(D5) 



(D6) 



D.2. SSC to synchrotron 

In this case we have to distinguish not only between the possible values of a, but also between the different 
cases of the Klein-Nishina-parameter K. 

For a ^ 1 we use for the synchrotron peak value equation (48 ), again, while for SSC peak value for iiT <C 1 
we use equation ( 57 ) . The ratio then becomes 

Rssc,.s = 4.4(1 -f l,,)a^. (D7) 



21 



For K ^ 1 we use equation ( 58 ) , giving 



Rssc.s — 26.8(1 + lec)-r^- 



(D8) 



Using equation ( |52[ ) for the synchrotron peak value in the case a 3> 1, we find with equation (601 the 
ratio in the cas 



Rssc.s = 6.8(1 + lec)a. 



The ratio for the case 1 ^ iiT <C a can be found by using equation ( 64 1 , yielding 

a 



5.2(1 + ?ec)7i 



74 



Lastly, we obtain the ratio for 1 a K with the help of equation (66): 

Rsscs = 14.4(1 + Zee) 7-^. 



(D9) 



(DIO) 



(Dll) 



Compared to ZS all ratios gain a factor (l + lec), raising the ratio, i.e. the Compton dominance, potentially 
by a lot. 

D.3. SSC to external 

Although this ratio is not between the two different components of the blazar SED, it still might be useful 
to give the ratios between the possible realizations of the high-energy component, since, depending on the 
parameters, it might be possible to discriminate between either the SSC or the external Compton scenario. 

D.3.1. Small injection parameter, a <C 1 

Firstly, we use the Thomson-limit for the external Compton SED, where the peak value is given by 
equation (70). Equation (57) gives the case K for the SSC component, yielding the ratio 

Rssc. 



1.6 a . 



For ^ 1 we obtain with equation ( 58 ) 



,(1 + U 



Rssc,ec — 14.3 ^ 3 ' 



(D12) 



(D13) 



For A' <C 1 we obtain with equation ( 57 ) : 



Secondly, we use equation (|7l|) for the peak value of the external Compton SED in the Klein-Nishina-limit. 

(D14) 



Rsscec = 6.3 X 10 ; « 7466 



For A' ^ 1 we use equation ( 58 ) , yielding 



R.sr.er. = 4.0 X 10 



5(1 -H/ec) a^te 



(D15) 



^Note the printing error in ZS. The correct value (apart from (1 + lec)) is given here. 
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D. 3. 2. Large injection parameter, a ^ 1 

In this case we obtain nine different ratios, which will be given below. 

Beginning with the case K <^lwe will use equation (601 for the SSC peak values. For v^r < v^^r we use 
equation ( 75 ) , which becomes 



Rs.c.ec = 1.0 X 10- 



(D16) 



In the case v^^r < Vbr < v-yor we take equation (77), yielding 



Rsscec = 8.0 X io2il±M (aj,e,,.f' . 



(D17) 



Thirdly, we have the case v^gr < where we use equation (81) to obtain 



Rsscec = ^A^^^a. 



(D18) 



Continuing with the case 1 <C if <C we take equation (64 1 for the SSC peak values. With equation 
(75) we get for v^r < v-ygr 



R - 7 5 X 10"*^^— t^-^ 

^ssc,ec — ( .O A lu ^ ^1/3' 



(D19) 



For v-yj^r < Vbr < i^7or we usc equation (77), again, giving 



R..r..n = 6.0 X 10- 



(D20) 



In the case i^jgr < i^br equation (81 ) gives the correct ratio: 



Rs 



6.3 



(1 + he) Q_ 

lec 



74 



(D21) 



Finally, we use equation (66) for the case 1 <C a <C i^. With equation (75) we yield for v^r < ^7 



Rssc.ec — 2.1 X 10 



(1 + lec) a^e-ec 



(D22) 



In the case i^-ygr < Vbr < v-^qt we take equation (77), obtaining 



R,,r..er. = 1-7 X 10" 



(1 + U a'/\li' 



I. 5/2 



(D23) 



Lastly, we have the case i^^^r < i^tr, where we find with equation (81 1 



R. 



17.9 



(1 + U 

he hi 



(D24) 



We see in all cases that for he ^ 1 the ratio is independent of this value, while for he ^ 1 the ratio becomes 
Rsse,ee oc 1^^- The latter implies a growing dominance of the SSC component over the EC component, which 
is expected. The former shows that it critically depends on the other parameters, which component might 
dominate the inverse Compton hump. 
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